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We discuss the evolution of anisotropic boost-invariant quark-gluon plasma possibly created at 
the early stages of relativistic heavy-ion collisions. Our considerations are based on the recently 
proposed formalism that is an extension of the relativistic perfect-fluid hydrodynamics. We analyze 
i) the pure partonic system described by the anisotropic phase-space distribution function and ii) 
the system of partons interacting with the local magnetic fields. The second analysis is a simplified 
attempt to include the effects of color fields on the particle dynamics. Our model results are 
discussed in the context of early thermalization and isotropization. Under general assumptions of 
the particle, energy, and momentum conservations we show that for large evolution times the ratio 
of the longitudinal and transverse pressures of the pure partonic system tends to zero. Hence, 
the system with the initial momentum distribution elongated along the beam axis always passes 
through the isotropic stage where the transverse and longitudinal pressures are equal. The inclusion 
of the magnetic field in this case gives negative contribution to the longitudinal pressure, hence the 
transient stage when the total longitudinal and transverse pressures become equal may be reached 
earlier, depending on the strength of the field. 
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I. INTRODUCTION 

The data collected in the heavy-ion experiments at the 
Relativistic Heavy-Ion Collider (RHIC) are most com- 
monly interpreted as the evidence that the matter pro- 
duced in relativistic heavy-ion collisions equilibrates very 
fast (presumably within a fraction of 1 fm/c) and its be- 
havior is very well described by the perfect-fluid hydro- 
dynamics [i|,M II i, i i [3 • 

Such features are naturally 
explained by the assumption that the produced matter is 
a strongly coupled quark-gluon plasma (sQGP) @. An- 
other explanation assumes that the plasma is weakly in- 
teracting, however the plasma instabilities lead to the fast 
isotropization of matter, which in turn helps to achieve 
equilibration Recently, it has been also shown that 
the model assuming thermalization of the transverse de- 
grees of freedom only is consistent with the data 
describing the transverse-momentum spectra and the el- 
liptic flow coefflcient V2 ■ This result indicates that the as- 
sumption of the fast thermalization/isotropization might 
be relaxed. 

In view of the problems related to the thermalization 
and isotropization of the plasma, it is useful to develop 
and analyze the models which can describe anisotropic 
systems. Recently, an effective model describing the 
anisotropic fluid/plasma dynamics in the early stages of 
relativistic heavy- ion collisions has been introduced 
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The model has the structure similar to the perfect-fluid 
relativistic hydrodynamics - the equations of motion fol- 
low from the conservation laws. However, it admits the 
possibility that the longitudinal and transverse pressures 
are different (as usual, the longitudinal direction is de- 
fined here by the beam axis). The main characteristic 
feature of the proposed model is the use of the pressure 
relaxation function R which determines the time changes 
of the ratio of the transverse and longitudinal pressures 
and, possibly, defines the way how the system becomes 
isotropic, i.e., how the transverse and longitudinal pres- 
sures become equilibrated. The role of the pressure re- 
laxation function is similar and complementary to the 
role played by the equation of state. It characterizes the 
material properties of the medium whose spacetime dy- 
namics is otherwise governed by the conservation laws. 

In this paper we develop the formulation of Ref. [llj 
in two ways: i) we introduce the microscopic interpre- 
tation for the relaxation function in the case where the 
considered system consists of particles whose behavior is 
described by the momentum-anisotropic phase-space dis- 
tribution function, ii) we include the effects of the fields 
considering the case of anisotropic magnetohydrodynam- 
ics - this approach may be regarded as a very crude at- 
tempt to include the effects of color fields on the particle 
dynamics. 

Our first finding is that the use of the anisotropic 
phase-space distribution functions leads inevitably to the 
pressure relaxation functions R which imply that the ra- 
tio of the longitudinal and transverse pressures tends 
asymptotically to zero, Pl/Pt 0. This behavior 
is complementary to the recent results obtained from 
the analyses of the early-stage partonic free-streaming 
[13, [ill- fa.ct, our approach includes the free-streaming 
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as the special case, however, it may be also applied in 
the cases where the collisions are present but their ef- 
fect does not change the assumed generic form of the 
phase-space distribution function. The time asymptotics 
Pl I Pt means that the systems with the initial pro- 
late momentum shape, i.e., the systems that are initially 
elongated along the beam axis in the momentum space 
with Pl > Pt (see for example Refs. [H, [M]), natu- 
rally pass through the transient isotropic stage where the 
transverse and longitudinal pressures are equal. Our sec- 
ond finding follows from the study of the magnetohydro- 
dynamic model. We show that the inclusion of the fields 
lowers the longitudinal pressure and increases the trans- 
verse pressure, hence, for the initially prolate systems the 
stage when the total longitudinal and transverse pressures 
become equal may be reached earlier, depending on the 
strength of the field. 

The two models discussed by us cannot explain the 
phenomenon of reaching the stable isotropic stage. How- 
ever, we indicate that the presence of the fields may have 
an impact on the process of isotropization, presumably 
restored by the effects of those collisions and/or field in- 
stabilities that are not taken into account in the present 
formalism. From the practical point of view, our for- 
malism allows for the determination of the space-time 
evolution of the color-neutral anisotropic distributions, 
which may be used, for example, as the background dis- 
tributions in the analysis of the plasma instabilities. In- 
terestingly, in the case where we have initially Pl > Pt, 
the dynamics of the background and the plasma insta- 
ble behavior evolve in the same direction, i.e., the two 
processes restore the equality of pressures. On the other 
hand, if the initial configuration has Pt > Pl, the plasma 
instabilities must compete with the growing asymmetry 
of the background (see Ref. [H, [l^ where the growth of 
the instabilities was studied in the Bjorken longitudinal 
expansion and substantially large times of about 20 fm/c 
were found for RHIC). The interplay of such competing 
processes may be directly related to the problem of very 
fast thermalization/isotropization taking place in rela- 
tivistic heavy- ion collisions. In addition, the discussed 
by us transformation of the longitudinal pressure into the 
transverse one is an interesting phenomenon analyzed in 
the context of the RHIC HBT puzzle - we note that the 
recently proposed explanations of this puzzle su gges t a 
very fast formation of the transverse flow fisl Il9j| . 

The paper is organized as follows: In Sect. II we 
consider the anisotropic system of partons described 
by the momentum-anisotropic phase-space distribution 
function. We calculate the moments of the distribution 
function, determine the pressure relaxation function R, 
and argue that the general form of R implies that the 
ratio of the longitudinal and transverse pressures tends 
asymptotically to zero. In Section HI we consider an 
example of boost-invariant magnetohydrodynamics. We 
analyze in detail the consistency of this approach and 
show that it may be treated as the special case of the 
formalism developed in pdj with the appropriate relax- 



ation function R. We also show how the presence of the 
local magnetic fields affects the dynamics of particles. 
We conclude in Sect. IV. 

Below we assume that particles are massless and we 
use the following definitions for rapidity and spacetime 
rapidity, 

1 gp+P|| 1 t + z 

y=oln^^ , ?7=-ln- , 1 

z i/p — p|| 2 t — z 

which come from the standard parameterization of the 
four-momentum and spacetime coordinate of a particle, 

= {Ep,p^,P\\) = {p±coshy,p^,pj_fimhy) , 
x'^ = {t,x±,z) — {TCOshr],x±,Tsinhri) . (2) 

Here the quantity p±^ is the transverse momentum 

Pl- = \JpI+pI, (3) 

and r is the (longitudinal) proper time 

T = - z2. (4) 

Throughout the paper we use the natural units where 
c = 1 and h = 1. 



II. ANISOTROPIC SYSTEM OF PARTICLES 

In this Section we consider a system of parti- 
cles/partons described by the distribution function which 
is asymmetric in the momentum space, i.e., its depen- 
dence on the longitudinal and transverse momentum is 
different. We calculate the particle current, the energy- 
momentum tensor, and the entropy current of such a 
system. As the special case we consider the exponen- 
tial Boltzmann-like distributions frequently used in other 
studies. This Section is also used to introduce general 
concepts of anisotropic plasma dynamics, which will be 
applied to the system consisting of particles and fields in 
the next Section. 



A. Anisotropic momentum distribution 

We take into account the phase space distribution func- 
tion whose dependence on the transverse and longitudi- 
nal momentum is determined by the two space-time de- 
pendent scales, and Ay, namely 

The form ([5]) is valid in the local rest frame of the plasma 
element. For boost- invariant systems, the explicitly co- 
variant form of the distribution function has the structure 
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where 



— {uq cosh rj^UxjUyjUQ sinh 77) , 



(7) 



1/^ = (smhr;,0,0,cosh7y), (8) 
and u^,Ux, Uy are the components of the four vector 

u^' ^ {u\u^,Q) ^ {u\ux,Uy,Q) . (9) 



The four-velocity is normahzed to unity 



1. 



(10) 



The four-vector 
plasma element. It 
Lorentz boost along 
pearance of the four- 
to the anisotropy - 
ement we have 
vectors and 
conditions: 



describes the four-velocity of the 
may be obtained from by the 
the z axis with rapidity 77. The ap- 
vector is a new feature related 
in the rest frame of the plasma el- 
= (0,0,0,1). We note that the four 
satisfy the following normalization 



1, V^^K 



-1, C/^K = 0. 



(11) 



The boost-invariant character of Eq. ([6]) is immedi- 
ately seen if we write the explicit expression for p ■ U and 
p ■ V which both depend only on y ^rj and the transverse 
coordinates, namely 

p-U p±uo cosh{y - rj) - pj^ ■ u±, 

p-V ^ p±uosmh{y - r]). (12) 



From Eqs. (|T2|) one can also infer that in the rest frame 
system of the plasma element, where 77 = and u± — 0, 
we have p ■ U — p± coshy and p ■ V — p± sinhy. Thus, 
in the local rest frame of the plasma element Eq. ([6|) is 
reduced to Eq. 



B. Moments of anisotropic distribution 

Using the standard definitions of and T^" as the 
first and the second moment of the distribution function 
([5]), namely 



d^p 



(27r)3 Ep 



(Pp 



(27r)3 Ep ' 

we obtain the following decompositions: 



(13) 



(14) 



(15) 



(16) 



We note that iV^ does not have the contribution propor- 
tional to the four-vector since such a term would be 
proportional to the scalar product V^N^ that vanishes 
in the local rest frame. Similarly, the energy-momentum 
tensor does not contain the terms proportional to the 
symmetric combination V^^U" + If-^V^, see Ref. [20| for 
a more explicit presentation of the analogous decompo- 
sitions. 

Equation (jlSp defines the particle density n, which may 
be calculated from the formula 



d^P J fp± bill 



(2^)3'' VAi' A|, 



(17) 





where we have introduced the dimensionless variables 



^11 = 



El 

All 



(18) 



In the similar way we calculate the energy density, 
d^p 



■Epf 



P± 



hj 

All 



(19) 



A2 A- 



27r2 





where the variable x is defined by the expression 

2 



X 



An 



(20) 



Finally, the transverse and longitudinal pressure is ob- 
tained from the equations 



Pt 



d^p 



Pa 



■f 



27r2 



(27r)3 2Ep 

dUCld^ 



PA_ hi 

A_l' All 



(21) 



/(a,cii: 



d^p P\\ fpi_ bill 
(2^)3 A|| 



(22) 



AiA^ rdUUdiwQ 



27r2 



-x& 



From now on the limits of the integrations over and 
^11 are always from to infinity. In the local rest-frame 
of the fluid element, where we have = (1,0,0,0) and 
V^' = (0,0,0,1) one flnds 



/ e 

Pt 

Pt 

VO Pl 



(23) 
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hence, as expected the structure ^ allows for different 
pressures in the longitudinal and transverse directions. 

One may also calculate the entropy current using the 
Boltzmann definition/ 



(27r)3 Ep \go 



1 - In ( ^ 

90 



(24) 



here go is the degeneracy factor related to internal quan- 
tum numbers such as spin or color. The entropy current 
has the structure 



(25) 



where 



Ah f 



1 - In 



50 



(26) 



Comparison of Eqs. (|17p and (j26p indicates that the par- 
ticle density and the entropy density are proportional, 
with the proportionality constant depending on the spe- 
cific choice of the parton distribution function /. 



It is quite interesting to observe that the structure of 
Eqs. (P?)) - (^5]) agrees with the structure derived in 
where no reference to the underlying microscopic picture 
was made but only the general consistency of the ap- 
proach based on the anisotropic energy-momentum ten- 
sor (jl6p and the conservation laws was studied. The only 
difference is that the entropy density a used in Ref. 
is now replaced by the particle density n. 

In fact, one may repeat the arguments presented in [Tlj 
replacing the assumption of the conservation of entropy 
by the assumption of the particle-number conservation 
(note that we have shown above that n and cr are pro- 
portional if one uses the ansatz In such a case we 
end up with the structure which exactly matches Eqs. 
(I?7)) - (pn)) and R may be identified with the pressure 
relaxation function. Moreover, the results of Ref. [ll| 
allow us to relate the variable x — A^/ A| with the quan- 
tity TiT^ - we shall come back to the discussion of this 
point below, after the analysis of some special cases of 
the anisotropic distribution functions. 



D. Boltzmann-like anisotropic distribution 



C. Pressure relaxation function 

With the help of the variables x — (A^/A||)^ and n we 
may rewrite our expressions (|19p . ((2T|) . and (|22p in the 
concise form 



Pt ^ - 



4/3 



4/3 



R{x) 



R{x) 



+ xR'{x) 



(27) 



(28) 



As the special case of the anisotropic distribution func- 
tion we may consider the exponential distribution of the 
form 



/i = 50 exp 



A 



Pi 

A2 



(32) 



'II 



which may be regarded as the generalization of the Boltz- 
mann equilibrium distribution where Aj_ = A|| = T (as 
explained above, go is the degeneracy factor connected 
with internal quantum numbers). In this case we recover 
the structure ([27| - (p9)) with the relaxation function of 
the form ^ 



Pl 



4/3 



R{x) 



2xR'{x) 



(29) 



where the function R(x) is defined by the integral 

_ -1/3 f d(.±^±d^\\ 



R{x) = X 



27r^ 



-^q+<lf{U.^\\), (30) 



Ri{x) 



3 50 X- 
27r2 



X arctan ^Jx — 1 



and the constant ([3T|) is simply 



50 

51 = :i2- 



(33) 



(34) 



R'{x) — dR{x)/dx, and g is a constant defined by the 
expression 



9= I ^^^/(^i'^ll)- 



(31) 



Another interesting anisotropic distribution function 
has the factorized form 



/2 = go exp ( exp f-y^ ) ■ (35) 



1 The formula I I24II assumes the classical Boltzmann statistics. It 
may be generalized to the case of bosons or fermions in the stan- 
dard way. 



^ Note that for x < 1 the function (arctan y/x — I) / y' x — 1 should 
be replaced by (arctanh \/l — x)/ ^1 — x 
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In this case we obtain 

-1/3 



2.0 



R2ix) 



27r2(l + j;)2 
3x 



1 + 5x^/x + 2x^\fx — 2x 



1 + VT+i 



and 



92 



go 

27r2' 



(36) 



(37) 



The calculation of the entropy density gives a — An and 
a — 8n, for the cases f — fi and f — f2, respectively. 

The structure of Eq. (f30|) implies that for x ^ I 
(A^ ^ A||) the function R{x) behaves like x^^^^. In 
this limit Pt — and e ^ Pl- Similarly, for a; 3> 1 
(A_L ^ A||) the function R{x) behaves like x^^^, implying 
that Pl = and e — 2Pt- This behavior is expected if 
we interpret the parameters A_l and A|| as the transverse 
and longitudinal temperatures, respectively. In agree- 
ment with those general properties we find 



R2{X) 



3go 
go 



-i/3 + i(in4-lna;)x2/3 
-^'^ + ^(61n2 - 8 - 31na;)x2/^ 



(38) 



for a; <C 1, and 



Ri{x) 

R2{x) 



3go 
47r 



,1/6 



50 { ^i/a 



-5/6 



-5/6 



(39) 



for x^ 1. 

In Fig. [T]we plot the ratios: Pl/Pt (solid line), Pl/s 
(decreasing dashed line), and Pt/s (increasing dashed 
line) for the two cases: / = /i (a) and f — f2 (b). 
The considered ratios are functions of the x parameter 
only. In agreement with the remarks given above we 
see that e = Pl for x = 0, and e = 2Pt in the limit 
X ^ oo. For f — fi the two pressures become equal 
if X = 1, since in this case the distribution function /i 
becomes exactly isotropic. For f — f2 the equality of 
pressures is reached for x « 0.7. Except for such small 
quantitative differences, the behavior of the pressures is 
very much similar in the two cases, as can be seen from 
the comparison of the upper and lower part of Fig. [TJ 



We note that the choice R{x) 



a/6 



corresponds to 



the case of transverse hydrodynamics, see Ref. In 
the transverse-hydrodynamic approach the matter forms 
non-interacting transverse clusters which do not interact 
with each other yielding Pl — 0. The concept of trans- 
verse hydrodynamics was initiated in Refs. ^U, [l^l and 
recently reformulated in Refs. [13, [13, [13 ■ 
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FIG. 1: (Color online) The ratios Pl/Pt (solid red lines), 
Pl/e (decreasing blue dashed lines), and Pt/e (increasing 
blue dashed lines) shown as functions of the variable x, a) 
the results for the distribution function (|32|) . b) the same for 
the distribution function (I35II. 



E. Time dependence of pressure anisotropy 

In this Section we briefly recall the arguments of Ref. 
[m concerning the consistency of the anisotropic plasma 
dynamics. The basic assumptions are the particle num- 
ber conservation, 

a^iV^ = d {nU") = 0, (40) 
and the energy-momentum conservation law, 

d^.T^" = 0, (41) 



with the energy- momentum tensor of the form (jl6[) . As 
shown in the previous Sections, the entropy conservation 
used in [ll| may be typically identified with the particle 
number conservation. In view of the further development 
of the model discussed in the next Section we shall turn 
to the particle number conservation as the basic input. 
We note that the assumption (1^0]) means that our de- 
scription may be valid only after the time when most of 
the particles is produced. 

The projection of the energy-momentum conservation 
law (|^T|l on the four-velocity Ui, indicates that the en- 
ergy density is generally a function of two variables. 
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e = e(n, r). The mathematical consistency of this ap- 
proach, i.e., the requirement that de is a total differen- 
tial, implies directly that the functions e{n, r), Prin, r), 
and PL{n, r) must be of the form (P7|) - ([25)) where 



X = Xq- 



(42) 



with xo, no and tq being constants that may be used to 
fix the initial conditions. In particular, it is convenient to 
regard tq as the initial time, and uq as the maximal initial 
density (at the very center of the system). Then, xq is the 
maximal initial value of x. Note that the particle density 
n is very small at the edge of the system, hence the initial 
transverse pressure is always zero in this region. On the 
other hand, at the center of the system at the initial time 
T = To we may have Pt < Pl or Pt > Pl depending on 
the value of xq. 

Combing Eqs. (PD)) and ([32]) we are coming to the main 
conclusion reached so far: For the microscopic phase- 
space distribution function of the form ^ , the pressure 
relaxation function is completely determined by Eq. (j30p 
where 



A- 



Xo 



nr 



(43) 



In the region where the matter is initially formed we have 
< n < TiQ and the right-hand-side of Eq. ((43|) grows 
with time - the particle density n cannot decrease faster 
than since this would require a three-dimensional 

expansion at the speed of light. We thus conclude that 
the ratio of the parameters A||/A_l tends asymptotically 
in time to zero. Consequently, for sufficiently large evo- 
lution times the ratio of the longitudinal and transverse 
pressures becomes negligible. As mentioned above, even 
if the initial conditions require that Pt is larger than 
Pl at the center of the system, at the edges we have very 
small density which means that P^ » Pt in this region. 
In the case where the longitudinal expansion dominates, 
n = nQTo/r and x — xqt'^/tq, hence x and r are simply 
related. 



F. Longitudinal free-streaming 

The anisotropic distribution functions considered in 
the previous Sections should satisfy the Boltzmann ki- 
netic equation (in some reasonable approximation). In 
this respect we assume that the effects of both the free- 
streaming and the parton collisions do not change the 
generic structure ([6]), while the time changes of the pa- 
rameters A^,A||, and are determined by the conserva- 
tion laws. The spirit of this approach is very similar to 
that used in the perfect-fluid hydrodynamics, where the 
collisions maintain the equilibrium shape of the distribu- 
tion function, whereas the conservation laws determine 
the time changes of the parameters such as temperature 
or the fluid velocity. 



Clearly the relation of our framework to the underly- 
ing kinetic theory should be elaborated in more detail 
in further investigations, which may possibly determine 
the microscopic conditions which validate our approxi- 
mations. Here, we may easily analyze the case of pure 
free-streaming where the distribution function satisfies 
the collisionless kinetic equation 



p^djix,p)^o. 



(44) 



For the pure longitudinal expansion (with vanishing 
transverse flow, u± = 0, and the parameters A_L,A|j de- 
pending only on the proper time r) we rewrite Eq. 
in the form 



cosh(y — T]) 



d sinh(y — rj) d 
dr T drj 



fiw,v)^0, (45) 



where w — pj_/X±{t) and v ^ p± sinh(j/ — 77)/A||(t), see 
Eqs. ([U and (fT2| . By direct differentiation we obtain ^ 

"dAii All 



df dXi_ ^ smh{y - t]) df 
An dv 



dw A^rfr 



dr 



= 0. 



(46) 



The solution to this equation exists for any form of the 
function / provided Aj_ is a constant and A|| ^ 1/r. Thus, 
we may write 




Tp± sinh(y — rj) 



where tq, A^ and Ajj are constants. 

In the considered case the variable x equals 



(47) 



A^ 
All 



(48) 



hence, it is consistent with Eq. (1431) . where for the 
boost-invariant longitudinal expansion we may substitute 
n = uoTo/t. We thus see that our approach includes the 
boost-invariant free-streaming as the special case. In par- 
ticular, Eq. (j47|) agrees with the form of the color neutral 
background used in Refs. [l^, [H, HJ- In the next Sec- 
tion we show that our framework includes also the case 
where the partons interact with local magnetic fields. 



III. LOCALLY ANISOTROPIC 
MAGNETOHYDRODYNAMICS 

In this Section we generalize the formulation discussed 
in Sect, [ill We analyze in detail magnetohydrodynamics 
as an example of the physical system consisting of par- 
ticles and fields, which is also known to exhibit strong 



^ For simplicity we consider here the functions depending on v'^ 
and disregard the sign of the absolute value. 
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anisotropic behavior. Of course, the magnetohydrody- 
namics by itself cannot be directly applied to modeling of 
the early stages of heavy-ion collisions. However, several 
phenomena analyzed in its framework show similarities 
with the color field dynamics discussed in the context of 
Color Glass Condensate [IE HH and Glasma [13] , hence 
we think that the elaboration of this example may shed 
light on more complicated color-hydrodynamics which 
may be the right description of the early stages of heavy- 
ion collisions. 

Our analysis of the boost-invariant magnetohydrody- 
namics, where the initial magnetic field is parallel to the 
collision axis, shows that in the considered system a sim- 
ilar phenomena take place as in the system consisting of 
particles only. The presence of the fields lowers the longi- 
tudinal pressure (which eventually may be negative) and 
increases the transverse pressure, see a related analysis 

in m. 



A. General formulation 

At first let us recapitulate the main physical assump- 
tions of locally anisotropic magnetohydrodynamics (for 
recent formulation see for example 29, 3^). Let be 
the plasma four- velocity and F'^'^ be the electromagnetic- 
field tensor. We define the rest-frame electric and mag- 
netic field by the following equations 



(49) 



(50) 



where is a completely antisymmetric tensor with 



^0123 



1. Eqs. gg]) and ^ yield 



We note that both i?^ and B^^ are spacelike and orthog- 
onal to [/^, 



E^E^ < 0, E^U^^O, 
B^.Bf' < 0, = 0. 



(52) 
(53) 



The picture of anisotropic magnetohydrodynamics re- 
quires that corresponds to the frame where the elec- 
tric field is absent. 



Ef" = 0. 



(54) 



In this case the Maxwell equations may be written in the 
form 



d*F^^'' = 0, 



(55) 



(56) 



where 

and is the dual electromagnetic tensor 



(57) 



(58) 



Besides Eqs. ([55]) - ([58l) the plasma dynamics is deter- 
mined by the particle conservation law, see Eq. (PD|) , the 
electromagnetic current conservation (the consequence 
of Eq. (|55p). and the energy- momentum conservation 
law for matter and fields. Before we analyze the role 
of the conservation laws we shall discuss, however, the 
constraints coming from the boost-invariance. 



B. Imposing longitudinal boost-invariance 

Our aim is to construct the boost-invariant field ten- 
sors F^^'^{x) and *F'^'^{x). The condition of boost- 
invariance requires that the transformed fields at the new 
spacetime positions are equal to the original fields in the 
new positions. Formally, this condition may be written 
in the form 



(59) 



where L describes the longitudinal Lorentz boost. Sim- 
ilarly, for a boost-invariant four- vector field A^{x) we 
have 



A^'ix') = L^'^A°'{x) = A'^ix'). 



(60) 



It is easy to check that the four-vectors and de- 
fined by Eqs. ^ and ([5]) are invariant under Lorentz 
boosts with rapidity a along the longitudinal axis, de- 
fined by the matrix 



I cosh a sinh a 

10 

1 

V sinh a cosh a 



(61) 



Since and are boost-invariant, the structure of 
Eqs. ((57|) and ((58)) suggests that the boost-invariant for- 
malism follows from the ansatz 



B^ = BV. 



(62) 



where i? is a scalar function depending on t and trans- 
verse coordinates x±. Equation (p^ defines the field 
tensors as the tensor products of the boost-invariant 
four-vectors, hence, by construction the field tensors are 
boost-invariant. In addition, we observe that in this case 
Eq. ([55)) is automatically fulfilled. 
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C. Homogeneous dual field equations 

Projection of the homogeneous dual field equations 
on the four-velocity Ui, gives 

Vd^.B + Bd^y - BUM^d^V = 0. (63) 

For the boost-invariant systems all terms in ([55)1 are iden- 
tically zero, hence it is automatically fulfilled. On the 
other hand, the projection of (I56|) on the four-vector 
gives 

C/^5^1n(^) -0, (64) 

hence B is related to the particle density n and the proper 
time r by the expression 

tl T 

B = Bo . (65) 

I 



The structure of the dual tensor allows us to infer the 
form of the electric and magnetic fields, 

BqUT I 1 Vi\ I \ 

B — (—UxSmhr],—UySimiri,u), (67) 

noTo 

E = (—""v cosh-q, Ux cosh?7, 0) , (68) 

The above structure implies directly that with no trans- 
verse expansion, i.e., for Ux = Uy = Q only the longitudi- 
nal magnetic field is present in the system and, in view of 
the relation n = uoTo/t, it should be a constant, B = Bq. 

In our general approach the situation Ux = Uy = cor- 
responds to the initial condition for the evolution of the 
system. It resembles the case of the Glasma [131 where 
also the longitudinal chromo-magnetic field is present, 
however, in the case of Glasma the direction of the field 
is random with the coherence transverse length set by the 
saturation scale (another difference is the presence of the 
longitudinal chromo-electric field in the Glasma). When 
the transverse expansion starts, due to the presence of 
the transverse pressure, it initiates the formation of the 
transverse magnetic and electric fields which are always 
perpendicular to each other, B ■ E = Q. We note, how- 
ever, that in the local rest frame of the plasma element, 
the only non-zero component is Bz- 

A more compact form representing the fields B and E 
may be achieved if we use the following parameterization 



One may check that with the ansatz ([65]) all four equa- 
tions in ([56|l are automatically satisfied, hence Eq. ((65|) 
is the main piece of information delivered by the homo- 
geneous dual field equations. In particular the equation 
V • B = turns out to be equivalent with the continuity 
equation for the particle number. 



Collecting now Eqs. ((62|) and ([65]) we find the explicit 
form of the dual field tensor 



—Ux cosh 77 I 
— UwCOshr; 

; 

(66) 

I 

of the particle current 

— n (u^ cosh ri,Ux,Uy,u'^ sinin]^ 

— (n u*^ coshyy, Tij;, n.y, rtM° sinhr;) . (69) 

Using the quantities Hx and Uy we write 



B = — — (—znx,—zny,nTu'^) , (70) 

"OTO 

E = -^{-tny,tnx,0). (71) 
noTQ 



D. Inhomogeneous field equations 

We turn now to the inhomogeneous field equations 
([55)1 . In our approach those equations may be used 
to determine the electromagnetic current of the sys- 
tem, j'' = {p, jxT jy, jz)- The straightforward calculation, 
where the form of the magnetic and electric fields given 



BqUT 

noTo 



Ux sinh T] Uy sinh rj 

-UxSmhr] 



-Uy sinh?7 











Ux cosh rj Uy cosh r/ 
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by Eqs. (I70p and ([71]) is used, leads us to the expressions 
Bot 



J =P = 



Jx 



- ]z 



noTQ 
Bp 

noTQ 
Bp 

npTp 

Bpz 

no To 



[rdyinu^) + 2ny + T^^-ny] , 
{dyUx - dxTiy) . 



(72) 



One may check by the expUcit calculation that the elec- 
tromagnetic four-current j'^ defined by Eq. (|72p is con- 
served, as required by the equation ([55]) . 

In the magnetohydrodynamic appoach one usually as- 
sumes that matter is neutral. In our case, the neutrality 
condition p — implies that the flow must be rotation- 
less, i.e., the following equation should be satisfied 

(73) 



dyU^ 



dxTly = 0. 



In this case also the longitudinal component of the elec- 
tromagnetic current vanishes, which means that the non- 
zero current circulates around the z-axis. 

The explicit calculation with the magnetic and electric 
fields given by Eqs. ([7(I| and ([71]) shows also that 



E + v X B :^0. 



(74) 



This is nothing else but the non-covariant version of the 
condition (|5T 



E. Conservation laws 
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FIG. 2: (Color online) The ratios: Pl/Pt (solid red lines), 
Pl/S (decreasing blue dashed lines), and Pr/i (increasing 
blue dashed lines) shown as functions of the variable x, 

a) the results for the distribution function ((32)) and c — 0.1, 

b) the same for c = 0.5. 



Besides the Maxwell equations, the equations of mag- 
netohydrodynamics include the conservation laws for: 
the particle number, the electromagnetic current (follow- 
ing directly from Eq. (|55p ). and the energy-momentum 
of the combined system consisting of matter and fields. 
The total energy-momentum conservation law may be 
written in the form 



(75) 



where the energy-momentum tensor, T'^'', including the 
contributions from matter and fields has the structure 



Pl-Pt 



An 
B^\ 



Pt 



(76) 



One may reduce the tensor ([76]) to the form ([T6|) if we 
introduce the following variables 

£ = e + — = e + e, 
Stt 



B 

Pt = Pt + ^=Pt 



Pl 



Pj. 



£2 
Stt 



Pj. 



Pl- 



(77) 



Clearly, the variables with a hat describe the sum of the 
matter and field contributions to the total energy density 
and transverse/longitudinal pressures (the field contribu- 
tions are marked with a bar) . 

Following the same method as that introduced in Ref. 
[Ill] we find that the energy-momentum conservation 
leads to the differential equation 



de 



Pt 



dn 



Pt-Pl 



dr. 



(78) 



Exactly this structure implies that the energy density 
and pressures are of the form (|27p - ([^5]) . So we may 
immediately write 



4/3 



R{x) 



4/3 



Pl 



\ 4/3 

9/ 



R{x) 



xR'{x) 



2xR'{x) 



(79) 



(80) 



(81) 
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where the complete relaxation function for matter and 
fields equals 

R{x) ^ R{x) +cax^^^, (82) 
with the parameter cq defined by the equation 



the equivalent form as 



r2 / „ \ 4/3 

^0 9 \ -2/3 

Co = — ' . 

Sn \noJ 



(83) 



In Fig. [5] we show the ratios: Pl/Pt (solid red lines), 
Ph/i (decreasing blue dashed lines), and Pt/s (increas- 
ing blue dashed lines) shown as functions of the variable 
X. Similarly to the case without the magnetic field we 
observe that the ratio of the longitudinal and transverse 
pressures decreases with x. The new feature of the case 
with the field is, however, that this ratio may become 
negative. This behavior reflects the negative contribu- 
tion of the field pressure Pl = —B'^/{8tt) to the total 
pressure Pl ■ It becomes dominant for the large values of 
X, where the matter contribution, growing as x^^^, may 
be neglected with the field contribution, growing as x^/^ . 

In view of our discussion in Sect. II, the variable x 
depends monotonically on time, hence the x dependence 
refiects to large extent the time evolution of the stud- 
ied ratios. If the initial conditions assume very small 
value of xq (and consequently of the initial x) the sys- 
tem has initially larger total longitudinal pressure than 
the transverse pressure Pt- The time evolution tends to 
equilibrate and then to invert the ratio of the two pres- 
sures. The time scale for this process is determined by 
the initial value of the field, Bq, as can be noticed by the 
comparison of the upper and lower part of Fig. [2l 

We close this section with the following remark. Since, 
_B is a function of n and r, we may rewrite Eq. (|78p in 



de 



Pl 



dn 



Pt - Pl 



B 



-dB. 



(84) 



This equation displays the dependence of the energy den- 
sity £ on the particle density n and the magnetic field B. 
The equation of the form e = e{n, B) plays a role of the 
equation of state. For the boost-invariant systems the 
functional dependence e(n, B) may be changed to the 
non-trivial dependence of e on n and r, as introduced in 
Ref. [13. 



IV. CONCLUSIONS 

In this paper we have developed the formalism intro- 
duced in Ref. [ll| discussing i) the system described by 
the anisotropic distribution function and ii) the system of 
partons interacting with local magnetic fields. The pre- 
sented results may be used to analyze anisotropic systems 
formed in relativistic heavy-ion collisions. In particular, 
they may be used to find anisotropic neutral distribu- 
tion functions which form the background for the field 
instabilities possibly responsible for the genuine thermal- 
ization/isotropization. In addition, our analysis indicates 
that the process of stable isotropization may require that 
the assumption concerning boost-invariance and/or en- 
tropy conservation should be relaxed. 
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St. Mrowczyhski for helpful discussions and critical com- 
ments. 



[1] P. F. Kolb, U. Heinz, in Quark- Gluon Plasma 3, edited by 
R.C. Hwa and X.-N. Wang (World Scientific, Singapore, 
2004), p. 634, |nucl-t h/0305084 
[2] P. Huovinen, in Quark-Gluon Plasma 3, edited by R.C. 
Hwa a nd X.-N. Wang (World Scientific, Singapore, 2004), 
p. 600, nucl-th/0305064 

[3] E. V. Shuryak, Nucl. Phys., A750 (2005) 64. 

[4] D. Teaney, J. Lauret, E. V. Shuryak, nucl-th/0110037 
[5] Y. Kama, et ai, Nucl. Phys., A774 (2006) 169. 
[6] C. Nonaka, S. A. Bass, Phys. Rev., C75 (2007) 014902. 
[7] T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey, Y. Nara, 

J. Phys., G34 (2007) S879. 
[8] E. Shuryak, J. Phys., G30 (2004) S1221. 
[9] S. Mrowczynski, Acta Phys. Polon., B37 (2006) 427. 
[10] A. Bialas, M. Chojnacki, W. Florkowski, Phys. Lett., 

B661 (2008) 325. 
[11] W. Florkowski, Phys. Lett., B668 (2008) 32. 
[12] W. Jas, S. Mrowczynski, Phys. Rev., C76 (2007) 044905. 
[13] W. Broniowski, W. Florkowski, M. Chojnacki, A. Kisiel, 

arXiv:0812.3393 [nucl-th]. 
[14] J. Randrup, S. Mrowczynski, Phys. Rev., C68 (2003) 



034909. 

[15] A. Rebhan, M. Strickland, M. Attems, Phys. Rev., D78 

(2008) 045023^ 

[16] A. Rebhan, arXiv:0810.2805' [hep-ph]. 

[17] W. Broniowski, M. Chojnacki, W. Florkowski, A. Kisiel, 

Phys. Rev. Lett., 101 (2008) 022301. 
[18] M. Gyulassy, Y. M. Sinyukov, I. Karpenko, A. V. 

Nazarenko, Braz. J. Phys., 37 (2007) 1031. 
[19] S. Pratt, arXiv:0812.4714 [nucl-th]. 

[20] R. Ryblewski, W. Florkowski, Phys. Rev., C77 (2008) 
064906. 

[21] U. W. Heinz, S. M. H. Wong, Phys. Rev., C66 (2002) 
014907. 

[22] U. W. Heinz, S. M. H. Wong, Nucl. Phys., A715 (2003) 
649. 

[23] M. Chojnacki, W. Florkowski, Acta Phys. Pol., B39 
(2008) 721. 

[24] M. Martinez, M. Strickland, Phys. Rev., C78 (2008) 
034917. 

[25] L. D. McLerran, R. Venugopalan, Phys. Rev., D49 
(1994) 2233. 



11 



[26] D. Kharzeev, E. Levin, M. Nardi, Phys. Rev., C71 (2005) 
054903. 

[27] T. Lappi, L. McLerran, Nucl. Phys., A772 (2006) 200. 
[28] J. Vredevoogd, S. Pratt, arXiv:0810.432 5 [nucl-th]. 



[29] M. Gedalin, I. Oiberman, Phys. Rev. E, 51 (5) (1995) 
4901. 

[30] M. Gedahn, Phys. Rev. E, 47 (6) (1993) 4354. 



